August 15, 2012 0:11 World Scientific Review Volume - 9.75in x 6.5in halo6 



3 



This chapter contains the following copyrighted materials: 

£N| . Reprinted Fig. 3 with permission from T. Nakamura et al., Phys. Rev. Lett. 96, 

o . 252502 (2006). C„ Pyrig M 2006 by the American Phvsica, Socie* 

£NJ , Reprinted Fig. 4 with permission from T. Aumann et al., Phys. Rev.C59, 1252 

kJTjf (1999). Copyright 1999 by the American Physical Society. 

Reprinted Fig. 1 with permission from K. Miernik et al., Phys. Rev. Lett. 99, 

192501 (2007). Copyright 2007 by the American Physical Society. 

Reprinted Fig. 2 with permission from A. Spyrou et al., Phys. Rev. Lett. 108, 
102501 (2012). Copyright 2012 by the American Physical Society. 

Reprinted Fig. 3 with permission from I. Tanihata et al., Phys. Rev. Lett. 100, 

192502 (2008). Copyright 2008 by the American Physical Society. 

Reprinted Fig . 1 with permission from A. Ozawa et al., Phys. Rev. Lett. 84, 5493 
£H | (2000). Copyright 2000 by the American Physical Society. 

Reprinted Fig . 2 with permission from I. Hamamoto, Phys. Rev. C69, 041306(R) 
(N 1 (2004). Copyright 2004 by the American Physical Society. 

>! 

(f) , Reprinted Fig . 2 with permission from P. Adrich et al., Phys. Rev. Lett. 95, 

132501 (2005). Copyright 2005 by the American Physical Society. 

Ln ■ 

Readers may view, browse, and/or download material for temporary copying pur- 
qq ■ poses only, provided these uses are for noncommercial personal purposes. Except as 

provided by law, this material may not be further reproduced, distributed, trans- 
mitted, modified, adapted, performed, displayed, published, or sold in whole or 
part, without prior written permission from the American Physical Society. 

>■ 

Reprinted Fig. 2 from Phys. Lett. B268, M. Fukuda et al., Neutron halo in Be 
studied via reaction cross sections, p. 339 (1991), with permission from Elsevier. 

Reprinted Fig. 3 from Nucl. Phys. A693, A. Ozawa, T. Suzuki, and I. Tanihata, 
Nuclear size and related topics, p. 32 (2001), with permission from Elsevier. 

Reprinted Fig. 2 from Phys. Lett. B707, M. Takechi et al., Interaction cross 
sections for Ne isotopes towards the island of inversion and halo structures of 29 Ne 
and 31 Ne, p. 357 (2012), with permission from Elsevier. 

Reprinted Fig. 3 from Phys. Rep. 264, T. Otsuka and N. Fukunishi, Nuclear mean 
field on and near the drip lines, p. 297 (1996), with permission from Elsevier. 



August 15, 2012 0:11 World Scientific Review Volume - 9.75in x 6.5in halo6 



2 



August 15, 2012 0:11 World Scientific Review Volume - 9.75in x 6.5in halo6 



Chapter 1 



Exotic nuclei far from the stability line 



K. Hagino 1 , I. Tanihata 2 ' 3 , and H. Sagawa 4 

1 Department of Physics, Tohoku University, Sendai 980-8578, Japan 
2 Research Center for Nuclear Physics (RCNP), Osaka University, Mihogaoka, 

Ibaraki, Osaka 5670047, Japan 
3 Research Center of Nuclear Science and Technology (RCN ST) and School of 
Physics and Nuclear Energy Engineering, Beihang University, No. 37 Xueyuan 
Road, Haidian District, Beijing, P. R. China 100191 
4 Center for Mathematics and Physics, University of Aizu, Aizu-Wakamatsu, 
Fukushima 965-8560, Japan 

The recent availability of radioactive beams has opened up a new era in nuclear 
physics. The interactions and structure of exotic nuclei close to the drip lines 
have been studied extensively world wide, and it has been revealed that unstable 
nuclei, having weakly bound nucleons, exhibit characteristic features such as a 
halo structure and a soft dipole excitation. We here review the developments 
of the physics of unstable nuclei in the past few decades. The topics discussed 
in this Chapter include the halo and skin structures, the Coulomb breakup, the 
dineutron correlation, the pair transfer reactions, the two-nucleon radioactivity, 
the appearance of new magic numbers, and the pygmy dipole resonances. 



1. Introduction 



Until the middle of 1980's, nuclear physics had been developed by investigating 
primarily stable nuclei which exist in nature. Many facets of atomic nuclei had 
been revealed, which include a mass, density distribution, radius, shell structure, 
collective excitations, and various decay modesP^ As a natural question, however, 
it had been discussed already in the late 1960's a question on how many neutrons 
can be bound for a given number of protonl^HH The first international symposium 
on nuclei far from the stability line was held in 1966 at Lysekil, Sweden,^ followed 
by the succeeding conference held in 1970 at CERNPThe questions which attracted 
nuclear physicists at that time include i) where is the neutron drip line? ii) do the 
nuclear models which were successful for stable nuclei work also for neutron-rich 
nuclei? and iii) the relevance to the r-process nucleosynthesis P It is remarkable that 
already in 1966 the 8 He nucleus was predicted to be stable by about 10 MeV against 
a dissociation to 6 He + 2n. 9 The n Li nucleus was also predicted to be slightly 
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Fig. 1. The matter radii of Li isotopes deduced from the measured interaction cross sections with 
a carbon target at 790 MeV/nucleon. The solid line shows the systcmatics known for stable nuclei, 
R oc A 1 / 3 , where A is the mass number. The experimental data are taken from Ref. 1101 

unbound, that is, only by 0.6 MeV, with respect to the 9 Li+2n configuration.^ 

The real start of the field of neutron-rich nuclei, however, was much later, that 
is, only in 1985, when an anomalously large matter radius of the n Li nucleus was 
experimentally discovered through the measurement of the interaction cross sec- 
tion.^ The matter radius was found to deviate largely from the known systematics 
in stable nuclei, which scales as A 1 / 3 as a function of mass number A of a nucleus 
(see Fig. 1). Together with a fact that the two- neutron separation energy, S^nj is 
extremely small for this nucleus (S2n was known to be 190± 110 keV at the time of 
1985p21 which has been updated to 378±5 keV^ 3 ), the large matter radius has been 
interpreted to be due to a long tail of the wave function for the weakly bound valence 
nucleonsP^ This structure is referred to as halo, in which the density distribution 
of valence neutron(s) largely extends over the core nucleus. Since the proton and 
neutron density distributions are almost the same in stable nuclei, the discovery of 
the halo structure was a big surprise. Subsequently, the interpretation of the halo 
structure was supported also by the observed narrow momentum distribution of the 
9 Li nucleus due to the breakup of n LiP2] It has been recognized by now that this 
exotic structure is an important characteristic feature of neutron-rich nuclei, and it 
has attracted much attention. 

The physics of neutron-rich nuclei has now been one of the main current subjects 
of nuclear physics. In fact, new generation RI beam facilities (such as RIBF at 
RIKEN in JapanJiS FAIR at GSI in Germany,^ SPIRAL2 at GANIL in France, 113 
and FRIB at MSU in the USA^) have been, or will soon be, in operation in the 
world wide. In this chapter, we summarize the developments of the physics of 
unstable nuclei in the past few decades. In Sec. II, we first discuss properties of 
one-neutron halo nuclei, that is, nuclei with a tightly bound core nucleus and a 
loosely bound valence neutron. We particularly discuss the role of single-particle 
angular momentum and the Coulomb dissociation. In Sec. Ill, we treat two-neutron 
halo nuclei. The main focus will be put on the correlation between the valence 
neutrons. As possible probes for the correlation, we discuss the Coulomb breakup, 




August 15, 2012 0:11 



World Scientific Review Volume - 9.75in x 6.5in 



halo6 



Exotic nuclei far from the stability line 5 

the two-nucleon radioactivity, and the two-neutron transfer reactions. In Sec. IV, 
we consider heavier neutron-rich nuclei. The topics to be discussed include the 
matter radii and the neutron skin thickness, the odd-even staggering of interaction 
cross sections, alpha clustering, the shell evolution and the deformation, and the 
collective excitations such as the pygmy dipole resonances. Finally, we summarize 
the chapter in Sec. V. 

As the physics of neutron-rich nuclei is diverse, it is almost impossible to cover 
all the topics in this chapter. We would like the readers to refer also to review 
articles EsHMl anc [ references therein. 

2. One-neutron halo nuclei 

2.1. Role of single-particle angular momentum 
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Fig. 2. The density distribution of the Be nucleus which is consistent with the measured inter- 
action cross sections (see the hatched area). Taken from Ref. 1361 



We first discuss properties of one-neutron halo nuclei, for which a weakly-bound 
valence neutron moves around a core nucleus. A typical example is n Be, for which 
the interaction cross section has been found considerably largep^S similar to the 
n Li nucleus. The one neutron separation energy, S n , is as small as 504± 6 keVp^ 
which can be compared to the one neutron separation energy of e.g., 13 C nucleus, 
$ n =4.95 MeV. This suggests that the 11 Be nucleus takes a halo structure. In fact, 
the large interaction cross section is consistent with a density distribution with a 
long tail, as has been shown in Ref. [36] (see Fig. 2). 

It is instructive to consider a simple two-body model with a core nucleus plus 
a valence neutron in order to understand the halo phenomenon. That is, we solve 
the Schrodinger equation, 



2m dr 2 



\)h 2 



2mr 2 



+ V(r) - e. 



ui{r) =0, 



(1) 
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Fig. 3. The radial wave functions u(r) for the relative motion between the valence neutron and 
the core nucleus 10 Be in the 11 Be nucleus. The upper panel shows the wave functions on the linear 
scale, while the lower panel shows the square of the wave functions on the logarithmic scale. The 
solid, the dashed, and the dotted lines correspond to the wave functions for the 2s, lp, and Id 
states, respectively. A Woods-Saxon shape is assumed for the mean-field potential, whose depth 
is adjusted for each angular momentum so that the single-particle energy is e = —0.5 MeV 
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Fig. 4. The root-mean-square (rms) radii for the 2s, lp, and Id states in the n Be nucleus as a 
function of the single-particle energy e. The meaning of each line is the same as in Fig. 3. The 
right panel is the same as the left panel, but plotted in the logarithmic scale for the horizontal 
axis. 



where m is the nucleon mass, I is the single-particle angular momentum, and V(r) 
is a spherical single-particle potential. Fig. 3 shows the radial wave function ui(r) 
on the linear (the upper panel) and the logarithmic (the lower panel) scales for the 
n Be nucleus with the two-body 10 Be+n model. We use a Woods-Saxon potential for 
V(r) with the radius and the diffuseness parameters of R=2.7A fm and a— 0.75 fm, 
respectively, whereas the depth parameter is adjusted for each angular momentum 
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I so that the single-particle energy e is —0.5 MeV. For simplicity, we do not consider 
a spin-orbit interaction. The solid, the dashed, and the dotted lines show the wave 
functions for the 2s, lp, and Id states, respectively. One can clearly see that the 
wave function for the 2s state is largely extended, while that for the Id state is 
spatially rather compact. The root-mean-square (rms) radii are 7.17, 5.17, and 
4.15 fm for the 2s, lp, and Id states, respectively. Fig. 4 shows the rms radii as 
a function of the single-particle energy. As one can show analytically,^ the rms 
radius diverges for I = and 1 (it behaves as |e| -1 ^ 2 for I = and |e| -1 ^ 4 for / = 1), 
while it converges to a constant value for higher values of I in the limit of zero 
binding energy. The halo structure, therefore, has been ascribed to an occupation 
of a weakly-bound I = or I = 1 orbit by a valence nucleon near the threshold! 37 * 38 ! 

2.2. Coulomb dissociation 

The halo structure significantly affects the dissociation process of a one-neutron 
halo nucleus in an external Coulomb field. The photoabsorption cross section for a 
dipole photon is given by 

_16vr 3 dB(El) 



where Ey is the photon energy and 
dB (El) 1 



dBL 2ji 



I KV/||eEiryi||Vi)r*(e/-e i -£l y ) ) (3) 



is the reduced transition probability (see e.g., Appendix B of Ref. [3]). Here, tpi 
and ipf denote the wave functions for the initial and the final states, respectively, ji 
being the angular momentum for the initial state ipi. eEi is the El effective charge, 
which is given by 

Z 1 A 2 -Z 2 A 1 
6E1= A 1 + A 2 e > (4) 

for a two-body system with a (A\, Z\) + {A 2 , Z 2 ) configuration. 

A characteristic feature of the dipole excitation is that the B(E1) distribution, 
dB(El) / dE, has a strong peak in the low energy region when the binding energy is 
small (that is, the soft El mode^tM). This large concentration of the El strength 
near the continuum threshold is caused by the optimal matching of wave functions 
between a weakly bound and continuum states. For a transition from an s-wave 
state to a p-wave state, the B(E1) distribution can be evaluated analytically as 

,3/2 



rfg(El) _ 2 y/lLp^ 

dE ~ rr 2 /! 6151 (\E b \+E c )^ [b) 
if one employs a Yukawa function for ipi and a spherical Bessel function for the 
radial part of ^ ^USHU] ^ gee R e f s [4"4l and HS1 for a general expression with arbitrary 
initial and final angular momenta). Here, [i is the reduced mass between the two 



August 15, 2012 0:11 



World Scientific Review Volume - 9.75in x 6.5in 



halo6 



8 



K. Hagino, I. Tanihata, and H. Sagawa 








■0.5 MeV 



1 MeV 







2 



4 6 
E c (MeV) 



S 



10 



Fig. 5. The B(E1) distribution from a bound s-wave state to continuum p-wave states given by 
Eq. (0 for two different values of the binding energy, E^. It is plotted as a function of the relative 
energy of the final continuum state, E c . 

fragments, Eb is the energy of the initial bound state, and E c is the energy of the 
final bound state, that is, the photon energy is E n — E c — Eb = E c + \Eb\. This 
equation indicates that the peak energy appears at E c = 3|Efc|/5 with the height of 
the peak being proportional to 1/|_E{,| 2 . The total strength is given by 



Therefore, the peak position moves towards low energy as the binding energy, \Eb\, 
decreases, and at the same time the height of the peak increases, leading also to 
the increase of the total El strength. These features can be clearly seen in Fig. 
5, which shows the -B(El) distribution given by Eq. ([5]) for two different binding 
energies. 

Notice that using Eq. (J3J) it is easy to derive that the total El strength (that 
is, the non-energy weighted sum rule) is proportional to the expectation value of r 2 
with respect to the ground state, 



As we mention in the previous subsection, the rms radius diverges in the zero 
binding limit for Z=0 and 1 states, leading therefore to a divergence of the total El 
strength. Thus, an observation of a large El strength makes a clear indication of a 
halo structure of the nucleus. 

Experimentally, the Coulomb dissociation of halo nuclei has been studied by 
Coulomb excitation experiments with a heavy target nucleus, such as 208 PbPSl The 
Coulomb breakup cross sections are often analyzed by the virtual photon theory, in 
which the cross sections are given as a product of the photo absorption cross sections, 
Eq. ©, and the virtual photon flux, N E1 (E 7 )^^ Fi gure 6 shows the experimental 
£?(E1) distribution for the 11 Be nucleus obtained by the Coulomb breakup reaction 
at 72 MeV/nucleonPs! The observed dipole strength shows a strong peak at about 
800 kcV excitation energy, which is consistent with the binding energy of about 





S(E1) = ^-e 2 E1 (r 2 )i. 



(7) 
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Fig. 2.6 in Ref. [28] 



Fig. 6. The experimental B(E1) distribution for the n Be nucleus deduced from the Coulomb 
breakup with a 208 Pb target at 70 MeV/nucleonPH Taken from Ref. [28] 

500 keV. The peak is large, and is again consistent with the halo structure of this 
nucleus. 

In addition to the large interaction cross sections} 49 * 50 ^ large Coulomb breakup 
cross sections are experimentally observed also for 19 C and 3 1 Ne .H3312] These nuclei 
are thus considered to be good candidates for halo nuclei. 

3. Two-neutron halo nuclei 
3.1. Two-nucleon correlation 

Let us now discuss properties of two-neutron halo nuclei, in which two valence 
neutrons are weakly bound to a core nucleus. For these nuclei, one must consider a 
(pairing) interaction between the valence neutrons. It has been well recognized that 
the pairing correlation plays an important role in nuclear physics JHHHH1 it leads 
to an extra binding for even-mass nuclei, and at the same time reduces the level 
density in the low energy region. Also, a pairing interaction scatters nucleon pairs 
from a single-particle level below the Fermi surface to those above, and consequently 
each single-particle level is occupied only partially. For weakly bound nuclei, the 
pairing interaction works by scattering nucleon pairs inevitably to unbound states. 

If one adopts a three-body model for a two- neutron halo nucleus, one could 
view it as a system of two interacting Fermions inside a confining potential. This 
is in a sense similar to a problem of interacting Fermion gas in a harmonic trap in 
atomic physics (see e.g., Refs. l56l and l57l and references therein). But a problem 
is much more challenging in weakly bound nuclei, because a trapping potential 
is not an infinite well (in contrast to a harmonic trap) so that the couplings to 
continuum arc important, and also because a trapping potential itself is constructed 
self-consistently from the interaction among nucleons. 
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Fig. 7. The two-particle densities for Li (the left panel) and for "He (the right panel) obtained 
with a three-body model calculation!^] These are plotted as a function of neutron-core distance, 
ri = r2 = r, and the opening angle between the valence neutrons, 812 ■ The densities are weighted 
with a factor 8ir 2 r 4 sin 612. 

Among two-neutron halo nuclei, the so called Borromean nuclei have attracted 
lots of attention. These are unique three-body bound systems, in which any two- 
body subsystem is not boundP^' Typical examples include 11 Li and 6 He, which 
can be viewed as three-body systems consisting of a core nucleus and two valence 
neutrons. Since both the n-n and n-core two-body subsystems are not bound, these 
nuclei are bound only as three-body systems. 

One of the most important current open questions concerning the Borromean 
nuclei is to clarify the characteristic nature of correlations between the two valence 
neutrons, which do not form a bound state in the vacuum. For instance, a spa- 
tial structure of two valence neutrons in the Borromean nuclei has attracted much 
attention. As a matter of fact, this has a long history of research as a general 
problem in nuclear physics. One of the oldest publications on this problem is by 
Bertsch, Broglia, and Riedel, who solved a shell model for 210 Pb and showed that 
the two valence neutrons are strongly clusterizedP^ Subsequently, Migdal argued 
that two neutrons may be bound in a nucleus even though they are not bound in 
the vacuum.^ The strong localization of two neutrons inside a nucleus has been 
referred to as dineutron correlation. It has been shown in Ref. 62 that the dineu- 
tron correlation is caused by admixtures of a few single-particle orbits with opposite 
parity. 

Although the dineutron correlation exists even in stable nuclei JsHMH it is en- 
hanced in weakly bound nuclei because the admixtures of single-particle orbits 
with different parities are easier due to the couplings to the continuum spectra. 
Probably it was Hansen and Jonson who exploited the idea of dineutron correlation 
explicitly for exotic nuclei for the first time. They proposed the dineutron cluster 
model and successfully analyzed the matter radius of n LiP^ They also predicted a 
large Coulomb dissociation cross section of the n Li nucleus. In the 1990's, more mi- 
croscopic three-body model calculations for neutron-rich nuclei startedpsEH These 
three-body model calculations have revealed that a strong dineutron correlation, 
where the two valence neutrons take a spatially compact configuration, indeed ex- 
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Fig. 8. The experimental B(E1) distributions for the 11 hi (the left panel) and for the 6 He (the 
right panel) nuclei deduced from the Coulomb breakup measurements. Taken from Refs. 1791 and 



ists in weakly-bound Borromean nucleiP^ 58 * 59 * 66 * ^ It has been shown that the 
dineutron correlation exists also in heavier neutron-rich nucleP^iHH] as we [\ as m 
infinite neutron matterP^^ The diproton correlation, which is a counter part of 
the dineutron correlation, has also been shown to exist in the proton-rich Borromean 
nucleus, 17 Ne!^ 

Figs. 7(a) and 7(b) show the two-particle density obtained with three-body 
model calculations for n Li and 6 He, respectively. These are plotted as a function of 
the neutron-core distance, r\ = T2 = r and the opening angle between the valence 
neutrons, 6*12. A weight of 47rr 2 • 27rr 2 sin#i2 has been multiplied. See Ref. [671 for 
the details of the calculations. One can see that a large fraction of two-particle 
density is concentrated in the region with small opening angle $12 for both the 
nuclei. This is a clear manifestation of the strong dineutron correlation discussed 
in this subsection. 



3.2. Coulomb breakup 

The Coulomb breakup of two-neutron halo nuclei can be discussed in a similar 
manner as that of one-neutron halo nuclei discussed in Sec. 2.2. The only difference 
is that the El operator is now given by 



D„ = e m RY lfl (R), 



where 



R 



(8) 



(9) 



is the center of mass coordinate for the two valence neutrons, and the El effective 
charge is given by 



2Z r 



eEi 



A c + 2 

with A c and Z c being the mass and charge numbers for the core nucleus. 



(10) 
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Fig. 9. The dipole strength distribution, d 2 -B(El) / de\de2, of 11 Li (the left panel) and 6 He (the 
right panel) as a function of the energies of the two emitted neutrons relative to the core nucleus. 
It is plotted in units of e 2 fm 2 /MeV 2 . See Ref. 82 for the details of the calculations. 

The left and the right panels of Fig. 8 show the measured B(E1) distribution 
for the 11 Li and 6 He nuclei E2M1 respectively. Those B(E1) distributions, especially 
that for the n Li nucleus, show a strong concentration in the low excitation region, 
similar to the B(E1) distribution for the 11 Be nucleus shown in Fig. 6, reflecting 
the halo structure of these nuclei. Moreover, the experimental data for n Li are 
consistent with the theoretical calculation only when the interaction between the 
valence neutrons is taken into account, strongly suggesting the existence of the 
dineutron correlation in this nucleus (see also Ref. ISTj) . 

The left and the right panels of Fig. 9 show calculated dipole strength distribu- 
tions, d 2 B{Fj1) / de\de2, of 11 Li and 6 He, respectively, obtained with the three-body 
modeP^l together with the Green's function method for the continuum dipole re- 
sponse.^ Here, e\ (e2) is the relative energy between the first (second) neutron 
and the core nucleus. One immediately notices that the strength distribution is 
considerably different between 11 Li and 6 He, despite similar ground state density 
to each other (see Fig. 7). This difference has been shown to be due to the different 
resonance properties of the neutron-core interaction between the two nuclei^ (see 
also Ref. 



3.3. Charge radii of halo nuclei 

Until recently Na was the lightest element of which charge radii (or proton density 
distribution radii) of neutron-rich short lived isotopes have been measured. Com- 
bined with the measurements of matter radii, the development of neutron skins 
thus has been presented for Na isotopes.^ For neutron halo nuclei, interaction 
cross section and fragmentation measurements have provided a mean to determine 
the matter density distribution. However the proton density distribution could not 
be determined directly from such measurements. Although the proton distribu- 
tion has been considered no to extend out even when neutron halos are formed, 
based on the fact that narrow momentum distributions, that indicate long tails of 
distribution, are observed only for a neutron or two-neutron removal channels of 
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fragmentation, no direct determination of proton distribution radii (or charge radii) 
was possible until recently. 

Recent developments of ion traps now provide a means to determine charge 
radii of very light neutron rich nuclei including neutron halo nuclei, 6 He, 11 Li, and 
n Be. The charge radii are determined by the isotope shift measurements of atomic 
transitions. In a very light atom, the isotope shift 8v includes two terms, 

8v = Sums + Su fs , (11) 

The first term Svms is the mass shift that is proportional to the difference of the 
masses (A and A') of two isotopes, 

A — A' 

d " MS 00 AA' ' ( 12 ^ 

and <5fFs is the field shift which is proportional to the difference of the rms radii of 
nuclei, 

Svfs oc Z x A[*(0)] 2 x 5(r 2 ), (13) 

where Z is the atomic number of the isotopes and A[5'(0)] 2 is the difference of 
the electron wave function at the nuclei. To obtain the field shift, the mass shift 
has to be calculated theoretically and then subtracted from the total isotope shift 
determined from the measurement. 

For a light nucleus as an example, the mass shift term for 6 He and 4 He is more 
than 10 5 times larger than the shift expected from a change of radius. The mass shift 
term cannot be separately measured experimentally so that an accurate theoretical 
estimation of this term is necessary in order to determine the charge radii. The 
recent development of atomic theory of a few electron system has enabled to do 
such a calculations successfully up to three electron systems P*l 

The charge radii of 6 He and 8 He have been determined by the ANL group 
using ANL/ AT LAS a nd GANIL RIB facilities to be 2.054±0.014 fm and 1.93±0.03 
fm, respectively! 86 * 87 ^ Those experiments use a magneto-optical trap of atoms for 
precision laser spectroscopy. The charge radii of Li and Be isotopes are determined 
by a GSI group using TRIUMF/ISAC facility and ISOLDE RIB facilities PH2l I n 
these experiments, a Doppler-free two-photon transition was used for the lithium 
measurement jSSL 9 ° l 91 l and collinear laser spectroscopy with a frequency comb for the 
beryllium isotopesP*! 

The proton distribution rms radii, R p = «/ (fp), has been calculated from the 

charge radii, y/ (r 2 ch ), by, 

(r 2 ch ) = (r 2 p ) + (pl) + ^(pl) + ^, (14) 

where r p is the radius of point proton distribution of a nucleus, p p and p n are 
the charge radii of free proton and free neutron, (p 2 ) — 0.769 ± 0.012 and (p 2 ,) = 
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Fig. 10. The matter, proton, and neutron rms radii (R m , Rp, and R n ) for He, Li, and Be isotopes 
determined from the interaction/reaction cross section and the charge radius measurements. 



-0.1161±0.0022fm 2 j22land the last term is the so called Darwin-Foldy term (0.033 
fm 2 )p«l rp^g re i a ^j on between matter, proton, and neutron rms radii is written as, 



ARi 



ZRi + NRt 



(15) 



The rms radii of matter (or point nucleon distribution), i? m , have been determined 
by the interaction cross section and reaction cross section measurements.^ The 
matter (Rm), proton (R p ), and neutron (R n ) radii so determined are presented in 
Fig. 10 for He, Li, and Be isotopes. One can clearly see that R p is much smaller that 
that of neutron R n in neutron rich isotopes and in particular in halo nuclei. However 
the proton radius increases slightly when neutron halo is formed. This observation 
is consistent with the view of a core+decoupled halo neutron(s) structure of halo 
nuclei that has been widely used for modeling halo nuclei. Under this model, the 
core of a halo nucleus has the same proton distribution as the isolated core nucleus. 
Because of the large distance between the core and the halo neutron(s), the core 
moves around the center of mass of the halo nucleus and therefore the proton radius 
in a halo nucleus is larger than that of isolated core. 
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Fig. 11. The geometry of the three-body model for two-neutron halo nuclei. 



3.4. Geometry of two-neutron halo nuclei 

Under the assumption of core+ two-neutron for a two-neutron halo nucleus, the 
spatial correlation of halo neutrons can be studied. The geometry of the model for 
such a nucleus is shown in Fig. 11. 

The relation between nucleon, proton, and neutron mean-square (ms) radii for 
a nucleus i is given by a similar equation as Eq. (|15p . Therefore from the measure- 
ments of matter and charge radii of a halo nucleus and its core nucleus, all radii 



of the halo nucleus, (r, n ), (rz), and (r„), and the core, 



{r%), and (r 2 n ), are 



determined. One can also define the ms matter radii of the halo neutrons, (r 2 ,), and 
that of the core, (r 2 m ), which are related as, 



A{r 2 m ) = A c (r 2 cm ) + A h (r 2 h ), 



(16) 



where A, A c , and Ah are the mass numbers of the halo nucleus, the core nucleus, 
and the number of halo neutrons, respectively. 

Using the coordinate of the core nucleus, p c , relative to the center of mass of 
the halo nucleus, that is, the movement of the center of the core in the halo nucleus, 
the relation between the nucleon, proton, and neutron radii between the core and 
the halo nucleus are represented as, 



2 ) 

C p I 



(r; 

(tin) 
(rln) 



(r 2 p ) = (rl p ) + (p 2 c ), 
(r 2 sm ) + (pl), 
{rln) + (Pc). 



(17) 
(18) 
(19) 
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From those equations, one can obtain the movement of the core, (p 2 )- Then the ms 
radius of the halo distribution, (r?), is also determined from Eq. (|16[) . 

The distance between the center of mass of the halo nucleus and the center of 
mass of the two halo neutrons (r^n) can be calculated from p c , 

A 2 h (rl n )=Al(p 2 c ). (20) 



\Zv2n) + V (Pc)- Let us define the size of di-neutron as the distribution radius of the 
neutrons around the center of mass of halo neutrons, r n _„ = r n i — r n2 = 2rdi_ n , 
and thus r 2 l _ n = 4r^_ , that is, rji-n is the radius of the di-neutron forming the 
halo. The radius of the halo distribution in the halo nucleus and the dineutron 
radius is related as, 

(rl) = (rL) + (4 i . n ). (21) 

Using these (r^) and {r 2n }, one can determine the di-neutron ms radius and thus 
the ms separation distance (v„_ n ) of the two halo neutrons. Notice that combining 
Eqs. (fl6|). (p]l, p0]>. and pi]) yields 



(r ™ } ~ I^2 (r ' m) + A~t~2 {j^2 {r2c - 2n) + 2 (r "- n> ) ' (22) 

for A h = 2. 

One can also obtain the two-neutron cross term (r„i • r n 2) as 

(r n i-r n2 ) = \(Al{pl)-{rl_ n )), (23) 

from 

Al(p 2 c ) = ((r nl +r n2 ) 2 ), and (r 2 n _ n ) = ((r nl - r„ 2 ) 2 ), (24) 

where Ah = 2 has been used. 

The empirical values for those variables for 6 He and 11 Li, extracted from the 
experimental interaction cross sections and charge radii, are shown in Table 1 and 
2. Corresponding theoretical values are also shown in the tables for a few model 
calculations. The da ta for 6 H e and n Li show that the distance between the core and 
the halo neutrons, (r^_ 2n ), is almost equal to the distance between two neutrons, 

sj ( r n— nli indicating that two-neutrons are sitting close together and have strong di- 
neutron correlations. The opening angle between the valence neutrons with respect 
to the core nucleus can be calculated from the empirical values for J (r 2 _ 2n ) anc ^ 

^(r 2 n _ n ) to be (0 12 ) = 58.9 ± 12.6 and (9 12 ) = 54.2 ± 3.04 degrees for n Li and 
6 He, respectively. (We should remark here that it is misleading to say that two 
neutrons are mostly sitting with opening angles obtained in this way. Instead, the 
mean opening angle is an average of a smaller and a larger correlation angles in the 
density distribution shown in Fig. 7.) 
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Table 1. The root-mean-square (rms) radii of 6 He. yj (r^J, \] (Tp), and \J (r 2 ) are the point nucleon, proton, and neutron 
rms radii of the nucleus, respectively. \J (rj£) is the rms radius of the halo-neutron distribution. U (r| n ) is the rms distance 



from the center of mass of the nucleus to the center of mass of the two valence neutrons, while y ( r c-2n) ' s * ne rms distance 
from the 4 He core to the center of mass of the two valence neutrons. The first estimation for the latter quantity by Wang 
et al. was 3.71±0.07 fmfSH yj (r 2 _ n ) is the rms distance between the valence neutrons, while {r„i ■ r n 2) is the correlation 
of the two valence neutrons. All these radii are given in the unit of fm, except for (r n \ ■ r n 2), which is given in the unit of 
fm 2 . See text for explanation. GFMC, NCSM, and AMD denote Greens Function Monte Carlo, No-Core Shell Model, and 
Anti-symmetrizcd Molecular Dynamics, respectively. 





Experiment 


GFMC 81 


VargaM 


Esbenscn 95 


Funada ub 


Zhukov 59 


NCSM UY 


ameIMI 




2.43± 0.03 




2.46 






2.45 




2.23 




1.912± 0.018 




1.80 








1.89±0.04 


1.83 




2.65± 0.04 




2.67 








2.67±0.05 


2.40 




0.808± 0.047 




0.87 














3.37± 0.11 
















yVi„) 


2.52± 0.05 




3.42 














3.84± 0.06 


3.81±0.20 




3.63 


3.51 


3.54 








3.93± 0.25 






4.62 


4.55 


4.58 






(r n i ■ r„2) (fm 2 ) 


2.70± 0.97 






0.54 


0.292 


0.325 







Table 2. Same as Table I, but for lx Li. The first estimation of ^(r 2 _ 2n ) by Sanchez 

et al. was 6.2±0.3 fm.EH SHF and TOSM denote Skyrme Hartree-Fock and Tensor 
Optimized Shell Model, respectively. 





Experiment 


Esbenserl^ 


SHF au 


Zhukov™ 


TOSM iuu 




3.50± 0.09 




2.87 


3.39 


3.41 




2.37± 0.04 




2.28 




2.34 




3.84 ± 0.11 




3.06 




3.73 




1.48± 0.12 












6.1± 0.3 












5.0± 0.5 












6.2± 0.5 


5.12 


6.26 


5.69 






7.0± 1.7 


6.77 






7.33 


(r n i ■ r n2 ) (fm 2 ) 


2.70± 0.97 











In principle, a similar analysis can be made for 8 He. However it is known that 
the possible core 6 He is known not to be a good inert core and thus 6 He + 2n 
model is very poor. Instead, it is known from the fragmentation experiments^ that 
4 Hc + An is a good model of 8 He ground state. Following the similar procedure 
as above, one can determine the movement of the 4- valence neutron center-of-mass 
in 8 He, \J {r\ n ) to be 1.07±0.05 fm. This value is much smaller than that in 6 He, 
V ( r 2n)/2 = 1-62 ± 0.03 fm. It is an indication that the four neutrons in 8 He are 
distributed more uniformly than in 6 He. 
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Fig. 12. An example of trajectories of the two emitted protons from the two-proton radioactivity 
of 45 Fe recorded by a CCD camera. Taken from Ref. 11141 



Notice that the core-2n rms distance, yf (R?) = J (r 2 ._2 n ) ) for two-neutron halo 
nuclei can be estimated also from the Coulomb dissociation cross sections. For the 
transition operator given by Eq. ([8]), the expectation value of R 2 in the ground 
state can be estimated from the total B(E1) value as 



(see also Eq. (f7|)). Using the experimental matter radii to estimate the rms distance 



y ( r ri-n)' the opening angle between the valence neutrons with respect to the core 
nucleus has been extracted as^^l (0 12 ) = 65. 2 ±12.2 degrees for n Li and 74.5 ±12.1 
degrees for 6 He (see also Ref. I103[) . which agree well with the results of the three- 
body model calculation of Ref. The value of (#12) for 11 Li is consistent also 
with the one obtained in Table 1 and 2, although the value for 6 He is somewhat 
larger. 

3.5. Two-nucleon radioactivity 

Although the experimental observation of the strong low-lying dipole strength dis- 
tribution in the n Li nucleus (see Fig. 8) has provided an experimental signature 
of the existence of dineutron correlation in this nucleus, it is still an open question 
how to probe it directly. That is, in the Coulomb breakup process, the ground 
state wave function of a two-neutron halo nucleus is perturbed by the external elec- 
tromagnetic field of the target nucleus, and it may not be easy to disentangle the 
dineutron correlation in the ground state from that in the excited states. It would 
be desirable if one could find an observable which reflects only the properties of the 
ground state. 

The two-proton radioactivity, that is, a spontaneous emission of two valence 
protons, of proton-rich nucleP^liMHloZ] is expected to provide a good tool to probe 





(25) 
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Fig. 13. The energy and the angular distributions of the two emitted neutrons from the two- 
neutron radioactivity of 16 Be. Taken from Ref. 11161 

the di-proton correlation in the initial wave function. Nuclei beyond the proton 

drip line are unstable against proton emission, but, since a proton has to penetrate 

the Coulomb barrier, their lifetime is sufficiently long to study their spectroscopic 

properties. A single-proton radioactivity has been found in many odd-Z proton-rich 

nuclei, and has provided a powerful tool to study the spectroscopy of proton-rich 

nuclei beyond the proton-drip li n e) 29 l 104 l 108 H llSI When the single-proton emission 

is energetically forbidden, proton-rich nuclei beyond the proton drip line decay via 

emission of two protons. Even though this process had been predicted theoretically 

in 1960p^ its first experimental discovery was much later, only in 2002p2HSl f or 

the 45 Fe nucleus. Subsequently, the energy and the angular distributions of the two 

emitted protons were also measured! 1 ^ An impressive development in Ref. 11141 

was the use of a new type of a gaseous detector, in which images of ionizing particle 

trajectories can be optically recorded with a CCD camera. An example of recorded 

trajectories of the two emitted protons from the two-proton radioactivity of 54 Fe is 

shown in Fig. 12. This technique was used also for the two-proton radioactivity of 
48 Ni ril5] 

Very recently, a ground state two-neutron emission was discovered for the first 
time for w Be^ This is an analogous process of the two-proton radioactivity, 
corresponding to a penetration of two neutrons over a centrifugal barrier. For the 
16 Be nucleus, the one-neutron emission process is energetically forbidden, which 
makes 16 Be an ideal two-neutron emitter. It is remarkable that the observed energy 
and the angular distributions of the two emitted neutrons show a strong indication 
of the dineutron correlation in the ground state of 16 Be (see Fig. I3)l!i£l 



3.6. Two-neutron transfer reactions 

It has been recognized for a long time that two-neutron transfer reactions are sensi- 
tive to the pairing correlationpJ^EJo] T/he probability for the two-neutron transfer 
process is enhanced as compared to a naive expectation of sequential transfer pro- 
cess, that is, the square of one- neutron transfer probability) 121 * 12 ^ The enhancement 
of pair transfer probability has been attributed to the pairing effect, such as the en- 
hancement of pair strength function^ 123 * 124 * and the surface localization of a Cooper 
pairP^EHl The pair transfer reaction is thus considered to provide a promising way 
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Fig. 14. The experimental angular distributions for the 1 H( 11 Li, 9 Li) 3 H reaction at 3 
MeV/nucleon. Taken from Ref. IT321 



to probe the dineutron correlation discussed in Sec. 3.1. However, the reaction dy- 
namics is rather complicated and has not even been well established. For instance, 
it is only in a recent calculation that a theoretical calculation achieves a satisfactory 
agreement with the experimental dataP^H It would therefore be not surprising that 
the role of dineutron in the pair transfer reaction has not yet been fully clarified. 

One example is a relative importance of the one-step (the simultaneous pair 
transfer) process and the two-step (the sequential pair transfer) process. In heavy- 
ion pair transfer reactions of stable nuclei, both processes are known to play a 
roleP^GjU p or weakly-bound nuclei, most of the intermediate states for the two- 
step process are likely in the continuum spectra. It is still an open question how 
this fact, together with the Q-value matching conditioners] alters the dynamics of 
the pair transfer reaction of neutron-rich nucleiP^ 

On the other hand, the cross sections for the pair transfer reaction of the Bor- 
romean nuclei, 11 Li and 6 He, have been measured recentlyP^2hl37| -j-lie data for 
the 1 H( 11 Li, 9 Li) 3 H reaction at 3 MeV/nucleon indicate that the cross sections are 
indeed sensitive to the pair correlation in the ground state of n Li (see Fig. 14)P^2I 
That is, the experimental cross sections can be accounted for only when the s-wave 
component is mixed in the ground state of 11 Li by 30-50%. Another important 
finding in this measurement is that significant cross sections were observed for the 
pair transfer process to the first excited state of 9 Lif^H which has made a good 
support for the idea of phonon mediated pairing mechanism .H3SI 

The two-neutron transfer from 6 He was investigated using 4 He and p tar- 
gets at FLNR(Dubna) at 151 MeV. | 133 | 134 ] The angular distribution from the 
4 He( 6 He, 4 He) 6 He reaction shows dominant contribution from di-neutron configu- 
ration (see the dashed line in Fig. 15 (a)). While the wavefunction contains both 
the cigar and di-neutron correlation the amplitude of the latter is larger (Fig. 15 (b), 
see also Fig. 7). The elastic scattering 4 He( 6 He, 6 He) 4 He was also studied at center 
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di-neutron cigar 



Fig. 15. (a) The angular distribution of the two-neutron transfer in the 4 He( 6 He, 4 He) 6 He reac- 
tionP^U The dashed line shows the calculation considering the dineutron component of the wave 
function while the dotted line is that with the cigar component, (b) The correlation densities in 
the wave function of 6 He for the dineutron and the cigar components. 

of mass energies of 11.6 MeV and 15.9 MeV at the ARENAS facility at Louvain-La- 
NeuveP^ An interpretation of the data at both energies^ 133 * 135 ! showed that the rise 
of cross section at large angle were due to two-neutron transfer. The calculation 
was carried out based on the coupled reaction channels approactP^ including both 
one-step and sequential two-step transfer with realistic form factors. It was further 
founcP^J that the direct two- neutron transfer strongly dominates over the sequential 
transfer at the low energies where the minimum in the angular distributior J 135 ! 13 2l 
is due to direct 2n transfer. The sequential transfer becomes more sizable for the 
higher-energy data, though direct 2n transfer still dominates. In the analysis in 
Ref. 11391 6 He was modeled as core +2n bound state wavefunction where two p 3 / 2 
neutrons coupled to J=0 (S-wave) and 1 (P-wave). The direct two- neutron trans- 
fer was found to be mainly due to the contribution from the S-wave. The P-wave 
part of the wave function belongs to the cigar type configuration, which is found to 
be a smaller contribution. The main conclusion is thus the same as Ref. 11331 and 
the dominance of S-wave two-neutron cluster transfer shows strongly correlated two 
neutrons in 6 Hc. 

The p( 6 He, 4 He)i angular distribution^^! (Fig. 16) when compared with the 
6 Li(p, 3 He) 4 He, shows a slightly larger cross section for 2n transfer at the center-of- 
mass (cm) scattering angles near 60 degrees. This was interpreted to be a signature 
of a more disperse 2n wavefunction in 6 He compared to the compact d in 6 Li. In 
both reactions a direct 2n transfer was assumed to be the dominant process. The 
p( 6 He,i) reaction was also studied at 25 MeV/nucleon at GANILp^u with an aim 
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Fig. 16. The a ngular distributions for p( 6 He, 4 He)t (the filled circles) and p( 6 Li, 3 He) 4 He (the 
open circles) P2H 

at looking into possible t+t cluster existing in 6 He beside the interest in exploring 
the He+2n configuration. The data are closer in magnitude and shape to the 2n 
transfer than the transfer of tritons. Therefore, the 4 He+2n configuration of 6 He is 
again confirmed. 

In the experiment with a heavier target, that is, the 6 He+ 65 Cu reaction at 
£'iab=22.6 MeV, both the 1-neutron (In) and the 2-neutron (2n) transfer cross 
sections were measured P^ An interesting observation for this system is that the 
cross sections for the 2n transfer are much larger than those for the In transferrin 
For stable nuclei, usually the In transfer cross sections are much larger than the 2n 
transfer cross sections and the opposite observation for the 6 He+ 65 Cu system 
can be regarded as a characteristic feature of a Borromean nucleus. A similar 
tendency, although less clearly, has been observed also for the transfer reactions 
for the 8 He+ 197 Au system at a similar energy £22] This measurement also shows 
that the transfer cross sections for the 8 He projectile are considerably larger than 
those for the 6 He projectile at energies around the Coulomb barrier J2Sl while these 
nuclei behave similar to each other in the subbarrier fusion reactions Further 
theoretical studies are apparently necessary in order to understand the differences 
and the similarities of these Borromean nuclei, 6 He and 8 He, in several reaction 
processes at energies around the Coulomb barrier. 



4. Heavier neutron-rich nuclei 

4.1. Matter radii and neutron skin thickness 



As we mentioned in Sec. 1, interaction cross sections 07 are intimately related to 
the size of colliding nucleipEH The interaction cross section is defined as a cross 
section for the change of the proton number Z and/or the neutron number N of a 
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Fig. 17. The left panel: Neutron rms radii for Na isotopes (the open symbols). These are 
deduced from the measured interaction cross sections using the proton rms radii (the filled circles) 
obtained with the isotope shift measurements. Taken from Ref. [841 The right panel: the density 
distribution for 132 Sn obtained with the Skyrme-Hartree-Fock method with the SIII parameter 
set. 

projectile nucleus after the interaction with a target nucleus. Another important 
quantity is a reaction cross section ap, which is defined as the total cross section 
minus elastic cross section, that is, the sum of the interaction cross section and the 
inelastic scattering cross section. Interaction cross sections are easier to measure 
than reaction cross sections, while the opposite is the case for a theoretical evalua- 
tion. For neutron-rich nuclei, cross sections for inelastic scattering are expected to 
be negligibly small at high incident energies liiSHlill anc j the interaction cross sec- 
tions are almost the same as the reaction cross sections. Because of this, measured 
interaction cross sections have often been compared to calculated reaction cross 
sections. 

Reaction cross sections have often been analyzed by the Glauber the- 
Or y| 30 | i43 | i45 | i46 | j n t ne optical limit approximation to the Glauber theory, together 
with the zero range approximation to the nucleon-nucleon interaction, the reaction 
cross section is given bySSEHl 
/' °° 

<tr — 2?r / bdb 
Jo 

where b is the impact parameter and s — (x,y) is the plane perpendicular to z. 
aNN is the total NN cross section, and p^ z \s) is defined by p( z \s) — j dzp(r), 
pp and px being the projectile and the target densities, respectively. It has been 
known that the optical limit approximation overestimates reaction cross sections for 
weakly-bound nuclei E32hlS2l j n order to cure this problem, Al-Khalili and Tostevin 
developed a few-body treatment for the Glauber theoryP^l Abu- Ibrahim and Suzuki 
have also proposed another simple method which effectively takes into account the 
higher order corrections 

The left panel of Fig. 17 shows the neutron rms radii for Na isotopes (the open 



1 — exp 



-aNN 



sp P 



b) 



(26) 
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symbols deduced from the measured interaction cross sections together with the 
proton rms radii (the filled circles) obtained from the isotope shift measurements. 
The figure indicates that the neutron rms radii are significantly larger than the 
proton rms radii for neutron-rich Na isotopes. This suggests that the neutron 
density distributions are largely extended over the proton density distributions, 
despite that there is no clear separation between a core nucleus and valence neutrons 
as in halo nuclei discussed in Sec. 2 and Sec. 3 (see Refs ll54t4159l for the halo 
structure in heavy nuclei). This structure is referred to as skin structure LLS2HI62I 
It should be mentioned that the skin structure is not necessarily related to s- and 
p- wave single-particle orbits, in contrast to the halo structure That is, the 
skin structure can be realized even with higher angular momentum states. As an 
example of skin nucleus, the right panel of Fig. 17 shows the density distribution for 
132 Sn obtained with the Skyrme-Hartree-Fock method^ 3 -' with the SIII parameter 
set.™ 

As we have discussed in Sec. 3.3, in order to discuss the skin thickness of neutron- 
rich nuclei, one needs both the matter and charge radii. For Na isotopes shown in 
Fig. 17, the proton radii have been obtained with the isotope shift measurements. 
The isotope shift measurement is not always applicable, however. In that case, one 
may use the charge changing cross section to estimate the rms radii for the proton 
distribution! 165 " 169 ! This is the cross section for a change in the charge number Z, 
and is considered to be sensitive to the proton distribution! 170 * 1 ^ Alternatively, 
one may use the proton elastic scattering measurement .UZ2UZ2I a yet novel method 
to extract the information on the proton distribution of neutron-rich nuclei is to 
use the electron scattering. Using a self-containing RI target (SCRIT), an attempt 
has already been successfully commenced at RIBF at RIKEN, JapanP^ 4 * 175 * 

4.2. Odd-even staggering of interaction cross sections 

The experimental interaction cross sections for neutron-rich nuclei often show a 
large odd-even staggering (OES). That is, the cross section for an odd-mass nucleus 
is significantly larger than the cross sections for the neighboring even-mass nuclei. 
A typical example is the interaction cross sections for 30 > 31 > 32 Ne, measured recently 
by Takechi et al. (see the left panel of Fig. 18). 50 In Ref. [1761 the odd-even 
staggering was analyzed by introducing the staggering parameter defined as 

/ , A cri(A + l)-2cj I (A) + f j I (A-l) 
73 = (-) g ' ( 27 > 

where 07 (A) is the interaction cross section of a nucleus with mass number A. The 
right panel of Fig. 18 shows the experimental staggering parameter 73 for Ne iso- 
topes as a function of the neutron separation energy for the odd-mass nuclei. It 
is compared with the results of the Hartrcc-Fock-Bogoliubov (HFB) calculations 
that takes into account the pairing correlations in the mean-field approximation. 
One can clearly see that the staggering parameter 73 increases rapidly for small 
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Fig. 18. The left panel: The experimental interaction cross sections for Ne isotopes with a 
carbon target at around 240 MeV/nucleon. The inset shows rms radii of stable nuclei. Taken 
from Ref. 1501 The right panel: the odd-even staggering parameter defined by Eq. II27I I for Ne 
isotopes as a function of the one-neutron separation energy for the odd-mass nuclei.EZSl The solid 
line shows the result of the Hartree-Fock-Bogoliubov calculations. 



separation energies, and goes up to a large value reaching 73 ~ 80 mb. Also, the 
experimental staggering parameters agree well with the HFB calculations, suggest- 
ing that the pairing correlation plays an important role in the odd-even staggering 
of the interaction cross sections. 



4.3. Alpha cluster in neutron-rich nuclei 

The a cluster model has been successful in describing the structure of light N = Z 
nuclei S ? 7 I 178 I In neutron-rich nuclei, extra neutrons are surrounding tightly-bound 
alpha particles. A new theoretical framework, the antisymmetrized molecular dy- 
namics (AMD), has been developed by Horiuchi et aZ.p^ which has been suc- 
cessfully applied to neutron-rich nuclei GBShlSSI j n t n j s method, a many-body wave 
function is assumed to be a parity- and angular momentum projected Slater de- 
terminant with multi-centered Gaussian single-particle wave functions. Generator 
coordinate method (GCM) calculations with such Slater determinants have also 
been considered!^] In this method, the alpha cluster is not assumed a priori, in 
contrast to the conventional cluster model, but it can emerge as a result of the 
energy minimization. The AMD calculations show that the cluster structure is de- 
veloped and stabilized in some neutron-rich nucleip^ which can be well interpreted 
in terms of the molecular-orbital picture. ^ 84 * 185 ! 
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Fig. 19. The measured _B(E2) values for TV = 20 even-even nuclei. The solid and the dashed 
lines are the results of the shell model calculation with the model space of sd+pf shells and sd 
shells, respectively. Taken from Ref. 196 



4.4. Shell evolution: change of spherical magic numbers in neutron- 
rich nuclei 

The shell closures and the associated magic numbers are one of the most important 
concepts in nuclear physicsP^ For stable nuclei, these magic numbers correspond 
to 2, 8, 20, 28, 50, 82, and 126. That is, if either the neutron number or the proton 
number (or both) coincides with one of these numbers, that nucleus is particularly 
rigid and takes a spherical shape. In general, the first 2 + state has a relatively large 
excitation energy in nuclei at the shell closures. However, a /3-decay spectroscopy 
studies had revealed by the middle of 1980's that the first 2 + state of the neutron- 
rich 32 Mg nucleus, which has N — 20, lies at as small as 0.886 MeV, suggesting a 
deformed shape of this nucleus JHZhlSai These experimental data stimulated lots of 
theoretical studies,^ 90 * 191 ! in which the nuclei around 32 Mg have been referred to as 
the nuclei in the island of inversionP^ That is, the relative position between the 
deformed 2-particle-2-hole intruder state and the spherical O-particle-0-hole state 
is inverted in this region. In 1995, Motobayashi et al. carried out the Coulomb 
excitation measurement of 32 Mg to the first 2 + state, and extracted a large B(E2) 
valueP^ The large B(E2) value is consistent with the nuclear deformation (see 
Fig. 19) JiSMHH and it has been concluded that the conventional N = 20 magic 
number does not hold in neutron-rich nuclei (see also Refs. 11971 and H98). A similar 
disappearance of shell closure has been observed also for N = 8^^ and N = 28p22 
A striking finding was that N = 16 becomes a spherical magic number in 
neutron-rich nuclei, accompanied by the disappearance of the N — 20 magicity. 
Ozawa et al. pointed out this fact by investigating systematically the neutron sep- 
aration energy in the neutron-rich p-sd and the sd shell regions^^ (see also Refs. 
I202H205| . When the neutron separation energy is plotted as a function of neutron 
number, it suddenly decreases across the shell closure. Ozawa et al. showed that it 
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Fig. 20. Systematics of one neutron separation energy, S n , for neutron-rich nuclei. The upper 
(a) and the lower (b) panels correspond to even-Z and odd-Z nuclei, respectively. The data for 
nuclei with same N — Z are connected with the lines. For the upper (lower) panel, with decreasing 
order, the lines correspond to N — Z = 1, 3, 5, 7, and 9 (N — Z = 0, 2, 4, 6, and 8). A small panel 
on r.h.s. shows spectroscopic factors of 2si/2 orbit in N=15 isotones. Taken from Ref. I201I 

indeed happens at N = 16 for nuclei with N — Z > 6 (see Fig. 20) indicating 
that the spherical magic number is shifted from N = 20 to 16 in neutron-rich nu- 
clei. Notice that the new magic number N = 16 implies that the neutron-rich 24 
nucleus with Z = 8 and N = 16 is a double magic nucleusP 06 1 207 1 

It is an important question to ask what makes the spherical magic numbers 
change in neutron-rich nuclei. Single-particle energies in a Woods-Saxon potential 
already show a quasi-degeneracy of 2si/ 2 an d ld 5 / 2 states and that the energy gap at 
N = 16 develops in weakly-bound systems! 201 * 208 ^ More microscopically, Otsuka et 
al. have argued that the tensor force as well as the spin-spin force of (<7-ct)(t-t) type 
have a responsibility for the change of shell structure in neutron-rich nuclei.^- 9 * 210 ! 
For instance, the attractive spin dependent interaction between the proton d 5 / 2 
orbit and the neutron d 3 / 2 orbits leads to a down-shift of the neutron d^/ 2 state in 
stable nuclei such as fijSiie, making the conventional shell gap at N = 20. In the 
s Oi6 nucleus, on the other hand, the protons do not occupy the d§/ 2 orbit, and 
thus this attraction does not work in the neutron orbits. Consequently, the energy 
shift of the neutron d 3 / 2 state does not happen, and the energy gap appears at 
N = 16 P U9 | 21U | A similar effect may affect also the shell structure in heavier regions, 
e.g., the magic numbers for superheavy elements ISliJ 



4.5. Deformed halo nuclei 



In the vicinity of the stability line, it has been well known that many open-shell nu- 
clei are deformed in the ground state. The nuclear deformation generates the collec- 
tive rotational motion, which is characterized by a pronounced rotational spectrum 
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Fig. 21. The probability of each angular momentum component in a single-particle wave function 
in a deformed mean-field potential. It is plotted as a function of the corresponding single-particle 
energy, (Lrn- 

Taken from Ref. [2181 

as well as strongly enhanced quadrupole transition probabilities. As we mentioned 
in the previous subsection, the first evidence for nuclear deformation in neutron- 
rich nuclei was the observation of a low-lying state in the 32 Mg nucleus. Another 
well known example is the so called parity inversion phenomenon in 11 BeP^ If one 
used a a naive spherical shell model, one would expect that the valence neutron in 
11 Be occupies the lp\/2 state in the ground state and the first excited states can 
be constructed by promoting the valence neutron to the sd shell, that is, the first 
excited state would be a positive parity state. However, the observed ground state 
of the n Be is 1/2+ state and the first excited state is 1/2" state at 0.324 MeV^i 
This appears as if the 2si/ 2 and the lpi/2 states are inverted in energy. This par- 
ity inversion problem has been naturally explained by considering that 11 Be is a 
deformed nucleus jHIlt2I£l 

A single-particle motion in a deformed mean-field potential is well known as 
a Nilsson orbit. 2 3 As a deformed mean-field potential does not have rotational 
symmetry, the corresponding single-particle wave functions are obtained as a linear 
combination of several angular momentum components. Misu et al. have pointed 
out that the s-wave component, when it contributes, becomes dominant in a de- 
formed wave function as the separation energy decreases, and eventually it gives a 
100% contribution in the limit of zero separation (see Fig. 21) l 2 ^ They have also 
shown that the p-wave component becomes dominant for negative parity states, 
although it does not give a 100% contribution even in the zero binding limit l 2 ^ See 
also Refs. 12181 and 12191 for related publications, and Refs P 20 1 2 -^ for self-consistent 
mean-field calculations for deformed halo nuclei. 

As we discussed in Sec. 2.1, the halo structure has been attributed to an oc- 
cupation of s or p single-particle orbit. The s and p wave dominance phenomenon 
in a deformed single-particle state suggests that the nuclear deformation enhances 
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a chance for a halo formation in weakly bound nuclei. The first evidence for a de- 
formed halo was observed recently for the 31 Ne nucleus p or this nucleus, both 
the Coulomb breakup cross section and the interaction cross section are found to 
be large, which clearly suggest a halo structure of this nucleus. Although a naive 
shell model indicates an occupation of the I/7/2 state by the valence neutron, the 
halo structure can be easily understood by considering a deformation, by which the 
p-w&ve component can be largely coupled in the wave function)222H224j p or ms tance, 
Urata et al. have carried out a particle-rotor model calculation^^ for this nucleus 
and shown that I* = 3/2~ is a good candidate for the ground state of 31 Ne, in 
which the component of the p 3 / 2 state coupled to the ground state of 30 Ne is as 
large as 44.9%P! 

4.6. Collective excitations of neutron-rich nuclei 

It has been well known that there are variety of collective excitations in atomic 
nuclei jUMHHHl p or stable nuclei, these collective motions are approximately clas- 
sified either as isovector or isoscalar types, in which the proton and the neutron 
motions are out-of-phase or in-phase, respectively. In neutron-rich nuclei, on the 
other hand, the isovector and the isoscalar modes are coupled to each other due to 
a large asymmetry in proton and neutron numbers! 227 * 228 ^ That is, a collective state 
may have both the isoscalar and isovector characters. In the extreme case, a pure 
neutron mode, in which only neutrons contribute to the collective excitation, may 
arise in neutron-rich nuclei.^- 9 * 230 ! A candidate for such neutron mode has been 
experimentally observed in 16 cP2HIII] 

There have also been lots of theoretical developments in descriptions of collective 
excitations in neutron-rich nuclei. Theoretically, the random phase approximation 
(RPA) has provided a convenient and useful method to describe excited states of 
many-fermion systems!^! In this method, excited phonon states are described as 
a superposition of many 1-particle 1-hole states. For weakly bound nuclei, the 
continuum effects play an essential role due to a much lower threshold energy com- 
pared to stable nuclei. The continuum RPA method was first developed by Shlomo 
and Bertschp23 which was subsequently applied to self-consistent calculations of nu- 
clear giant resonances with Skyrme interaction by Liu and Van Giai! 2 ^ Hamamoto, 
Sagawa, and Zhang have extensively applied this method to nuclear responses in 
neutron-rich nucleiP^ 7 * 228 * 236 * 2 -^ An extension of the continuum RPA to deformed 
nuclei has also been carried out by Nakatsukasa and YabanapSHl Another important 
development was to include the pairing effects in the continuum RPA. As we dis- 
cussed in Sec. 3, the pairing and the continuum couplings play an essential role in 
neutron-rich nuclei. These effects can be taken into account by extending the RPA 
to the quasi-particle RPA (QRPA). The continuum QRPA method on top of the 
ground state described by the Hartree-Fock-Bogoliubov method has been developed 
by several groups 1I 23 I 124 I 230 I22SM24Z| anc j nag been applied to neutron-rich nuclei. 

Recently, much attention has been paid to low-lying dipole (El) strength in 
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Fig. 22. The experimental data for the Coulomb excitation (the left panels) and the photo- 
neutron (the right panels) cross sections for 124 Sn, 130 Sn, and 132 Sn. Taken from Ref. 12541 

neutron-rich nuclei, which has often been referred to as pygmy dipole resonance. It 
has been well known that the dipole strength is far dominated by the giant dipole 
response (GDRjPMHU ( see a j so p e f |248l for a recent complete measurement for 
the dipole strengths in 208 Pb over a wide range of excitation energy). However, 
nuclei with a neutron excess often show a low-lying dipole strength at energies 
much lower than the GDR. The pygmy dipole resonances have been experimentally 
found not only in neutron-rich nuclei^HHH b u t also in stable nucleP^SMsey ( see pjg_ 
22). The pygmy dipole resonance is important also from the astrophysical point of 
view, as it affects significantly a radiative neutron capture rate, which is relevant 
to the r-process nucleosynthesis!^ Although the exact nature of the pygmy dipole 
resonance has not yet been clarified completely, it has been pointed out that the 
pygmy dipole strength is strongly correlated with the neutron skin thicknes d 268 * 269 ! 
(see, however, also a counter argument in Ref. I270p . As the neutron skin thickness 
is intimately related to the equation of state (EOS) in asymmetric nuclear mat- 
ter p22 it is expected that the information on the nuclear matter properties, such as 
the symmetry energy coefficients, may be obtained by studying the pygmy dipole 
strength in neutron-rich nucleiP^l 

5. Summary 

Physics of unstable nuclei has been developed rapidly thanks to the recent availabil- 
ity of radioactive beams in the world. A new era has been commenced in nuclear 
physics, where the isospin degree of freedom, that is one of the fundamental quan- 
tum numbers in atomic nuclei, can be controlled in self-bound interacting Fermion 
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systems. Many new features of atomic nuclei have been discovered or theoretically 
discussed so far, as we have reviewed in this Chapter. We list some of these below: 

(1) The most prominent discovery was the halo structure, in which the neutron 
density distribution largely extends over the proton distribution. In stable 
nuclei, the neutron and the proton distributions arc similar to each other within 
a scaling factor, and thus the halo structure is a clear manifestation of the 
decoupling of proton and neutron in weakly bound nuclei. The halo structure 
is ascribed to an occupation of either s or p orbits, for which the root mean 
square radius diverges in the zero binding limit. Typical examples of the halo 
nuclei include 6 He, 11 Li, and 11 Be. 31 Ne has also been identified recently as a 
deformed halo nucleus. 

(2) A similar decoupling effect is the skin structure. While the halo structure 
corresponds to a long low-density tail of neutron distribution, the skin structure 
corresponds to a layer of extremely neutron-rich matter. The matter radii 
of neutron-rich nuclei have been systematically studied with interaction cross 
section and isotope shift measurements. 

(3) Neutron-rich nuclei often show a soft dipole excitation, that is, dipole strengths 
in the low excitation energy region. For halo nuclei, the soft dipole mode 
is attributed to the threshold effect, that is, the optimal matching of wave 
functions between a weakly bound and continuum states. For skin nuclei, the 
low-lying dipole resonances have been referred to as the pygmy dipole mode. 
Even though the exact nature of the pygmy dipole mode has not yet been fully 
clarified, it has been expected that it is closely related to the equation of state 
of asymmetric nuclear matter. 

(4) The dineutron correlation in neutron-rich nuclei, such as n Li and 6 He, has been 
theoretically predicted. This is a spatial correlation, with which two valence 
neutrons take a compact configuration. The recent experimental data for the 
Coulomb dissociation of 11 Li strongly suggests the existence of the dineutron 
correlation in n Li. A similar correlation in proton rich nuclei, that is, the 
diproton correlation has also been predicted. The two-nucleon radioactivities 
as well as the two-neutron transfer reactions are expected to provide a direct 
probe of the dineutron correlation. 

(5) Another important feature in neutron-rich nuclei is the change of (spherical) 
magic numbers. It was considered that the magic numbers exist independently 
for proton and neutron. However in neutron-rich nuclei some of the magic 
numbers have been found to disappear (N=8 and 20) and a new spherical 
magic number appears (A=16). These changes have revealed the importance 
of tensor interaction. 

A new generation RI beam facility, RIBF, at RIKEN, Japan, has already been in 
operation, and other new generation facilities, such as FAIR (Germany), SPIRAL2 
(GANIL), and FRIB (USA), will also be in operation in a few years. We are now 
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at a stage in which we gain a deep insight in nuclear many-body systems, from 
stable nuclei to weakly-bound unstable nuclei in a unified manner, and are about 
to explore a vast terra incognitcP^m the nuclear chart. 
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